Quark recombination and elliptic flow 



Scott PratlQ 

Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824 

Subrata PaQ 

National Superconducting Cyclotron Laboratory and Department of Physics and Astronomy, 
Michigan State University, East Lansing, Michigan 48824 

Elliptic flow systematics for different hadron species have been explained by quark coalescence 
models. It has been argued that the elliptic asymmetry v% should scale with the number of quarks 
that comprise the hadron. We show how these arguments are sensitive to the relative role of 
asymmetries in the phase space density vs. asymmetries in the effective volume of emission. We 
also discuss the degree to which coalescence arguments differ from thermal models. Illustrative 
calculations based on solving the Boltzmann equation are presented along with the results of blast- 
wave models. Although the issue is complicated, ambiguities might be clarified by measurements of 
source-size parameters for nucleons at higher transverse momenta. 
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One of the surprising early measurements from RHIC shows that the proton/pion ratio reaches or even exceeds unity 
for transverse momenta p± above 2 GeV/c One explanation for this phenomena is that quarks originating from 
different nucleon-nucleon collisions recombine via coalescence mechanisms & El If this is indeed the explanation, 
, it will make a case for having created a new state of matter where memory of the original color singlet hadrons is lost 
as the constituent partons become deconfined before going through a "mix-and-match" process, referred to as quark 
recombination. 

Coalescence arguments have also been invoked to explain systematics of elliptic flow asymmetries 0, 0, HI El 
These analyses are based on the parameter H| , 

3' , . fd<f>(dN/d<f>dp t ) cos 20 



where the azimuthal angle <f> is measured relative to the reaction plane. It is argued that if the underlying quarks 
have an asymmetry v^ 1 , hadrons will have an asymmetry 0- 0- 0] ■ 

v^{p t )=n q v^(p t /n q ), (2) 

where n q is the number of quarks comprising the hadron. This might explain the fact that baryons seem to have 
~ 50% higher asymmetries than mesons 0]. 

The main purpose of this paper is to investigate the validity of the scaling relation for V2 of Eq. Aside from the 
usual caveats concerning coalescence, Eq. requires an assumption that the average phase space density has the 
same asymmetry as the invariant spectra. After demonstrating how this assumption plays a role we present sample 
calculations to illustrate the circumstances under which Eq. @ is valid. For blast-wave models of the emitting 
state 0, ^2 ; we will show that the validity is sensitive to seemingly arcane choices of how to parameterize 

the model. Boltzmann equations are also solved for the partonic evolution where validity of Eq. (0 is found to be 
sensitive to the initial conditions. 

The next section presents a discussion of quark recombination through coalescence with an emphasis on explaining 
the subtle differences between coalescence and thermal models. The subsequent sections present results from the 
blast-wave and Boltzmann calculations respectively. Results are summarized in Sec. IVl 
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II. THEORY: COALESCENCE AND ELLIPTIC FLOW 

Coalescence arguments have been applied in numerous instances |3 EH El> 0> HE El> to problems of fragment 
production in nuclear physics. It has been especially successful in describing the systematics of light nuclei production 
in high energy heavy-ion collisions where a high degree of thermalization is attained and temperatures far exceed 
typical nuclear binding energies. Coalescence shares many features of thermal models, and in fact provides the same 
results in the limit of high temperature relative to the binding energy pl| . Although the focus of this paper is quark 
coalescence, the theory discussed in this section can be equally well applied to the coalescence of nucleons to form 
deuterons or light clusters. The main theoretical difference with quark coalescence is that one must be more cognizant 
of the effects of the non-zero binding energy and of relativistic motion. Whereas the deuteron binding energy is only 
2.2 MeV and the momentum components of the nucleon-nucleon relative wave function has most of its strength in 
the range of tens of MeV/c, constituent quarks may form hadrons whose mass differs from the sum of the constituent 
quark masses by many hundreds of MeV. Momentum components of the quark's relative motion also tend to be in the 
range of hundreds of MeV/c. Furthermore, one must add the qualifier that hadrons cannot necessarily be expressed 
in terms of constituent quarks, which ignores the gluons, sea quarks or vacuum fluctuations of meson fields. 

In this section we further outline some of the subtle issues mentioned above, but rather than focusing on these 
issues, we will mainly consider a simplified case where n q constituent quarks each of mass m q combine to form a 
hadron of mass n q m q . This allows us to sidestep some of the more delicate points about coalescence and derive some 
simple relations for the behavior of elliptic flow as a function of the number of quarks. 



A. Review of Coalescence Theory 

The meaning of the term coalescence varies somewhat throughout the literature 0, 0, 0, IT?! fill fill l2fij | . For the 
purposes of our discussion, we assume that it refers to the formation of a bound state of two pre-existing constituents 
via the sudden approximation, i.e., at the moment of coalescence the binding interaction is turned on suddenly and 
the probability that constituents a and b will form a composite object C is 

/c(Po,R,tc) = j d 3 rd 3 qW^\q 7 r)f a { ma Pc/{m a + m b )+q,R + r/2,t c )f b {m b P c /{rn a + m b )-q } Il-r/2,t c ). (3) 

Here, the phase space densities are denoted by /, and W refers to the Wigner decomposition of the relative squared 
wave function. Pp is the momentum of the composite object, t c the coalescence time, and R and r are the center-of- 
mass and relative coordinates. Although this expression is non-relativistic, the problem can be considered in a frame 
where the velocity of the constituent particle V is zero or small. 

If the characteristic spread of the relative momenta and relative coordinates in VF(q, r) are small compared to the 
scale of the initial momenta and spatial sizes, W may be replaced with a delta function, 

^ ) (q,r)«5(r)5(q). (4) 

The simplified coalescence formula is then, 

/c(Pc,R, tc) « f a (m a Pc/(m a + m&),R,i c ) f b (m b P c /(m a +m&),R,t c ). (5) 

For a thermal source, we consider a Boltzmann form for the phase space density, 

/i(miV,R,t c ) = exp [-(£.; - Mi (R,t c ))/T(R,t c )] , (6) 

where E 1 refers to the energy in the frame of the thermalizing medium. Inserting the thermal forms for f a and f b in 
Eq. (J5J nearly reproduces the thermal form for fc with fic = Ha + Hb- It differs by the factor exp(—B'/T) where B' 
is the binding energy measured in the frame of the thermalized medium, 

B' = E' a + E' b -E' c . (7) 

For the coalescence of nucleons into deuterons, the binding energy of 2.2 MeV is much smaller than the temperature, 
making the coalescence and thermal forms nearly indistinguishable. 

The source- function formalism, which is commonly used in two-particle correlations |22j | , can also be applied to this 
problem. In this formalism the last scatterings with third bodies are treated as randomizing interactions which can 
place the particles into bound states. The source functions are defined in terms of the quantum T matrices which 
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describe the scattering from third bodies into the final two-body state. The momentum distribution of a composite 
object created from two particles, 



d 3 p, 



•c 



E 



d 4 x a d 4 x b Tf,(x a ,x b )U(x a ,x b ;pc) 



(8) 



where U(x a , x b ;pc) is the evolution matrix for particles evolving from space-time points x a and x b into the asymptotic 
state of a bound particle, C, with momentum pc- The remainder of the system will evolve into the state /'. By 
imposing the condition that the sources are independent, (apart from the final-state interaction), the T matrices can 
be expressed as a product, 



T f ,(x a ,x b ) -> T fk (x a )TJ,(x b ), ~> E ' 

/' fUl 



(9) 



where f' a and f b refer to the two independent sources of the particles a and b. Inserting the factorized form into Eq. 



dN^ 
d 3 p c 



= E 

faJi 

= E 

/a, A 



d 4 x a d 4 x b T fa (x a )T fb (x b )U(x a ,x b ;pc) 



d 4 x a d 4 x b Tf a (x a )T fb (x b ) exp [i(p a + p b ) ■ (m a x a + m b x b )/ (m a + m b )} u c (x' a - x' b ) 



(10) 



where the two scattering centers evolve into states f' a and f b . The center-of-mass behavior has been factored out 
of the evolution matrix, and u is the evolution matrix in the frame where Pc = 0, with X' a ~ X' b referring to the 
coordinates measured in that frame. (If the sources are not independent, there would have been correlations in the 
emission of a and b in the absence of final-state interactions.) 
Invoking properties of Fourier transform, 



dN^ 



d 4 X a d 4 X b d 4 rd 4 q S a (m a Pc/ (m a + m b ) + q, X a )S b (m b Pc/(m a + m b ) - q, X b )e~ 



d 3 P c (2tt) 

x u*{X' a -X' b + r'/2)u(X' a -X' b - r'/2); 

Si( Pi ,X) = E/ cPxT^X + x/2)T fi (X - x/2)e- ipi - x . 



(11) 



When the source function Si(p, x) is evaluated on-shell, po — E p , it can be identified with the probability of emitting 
a particle of type i with momentum p from space-time point x. 

A more tractable expression can be obtained by applying the smoothness approximation p^ . l2^ | , 

S a (m a Pc/ (m a +m b )+q, X a )S b (m b P c / (m a +m b )-q, X b ) sa S a {m a P c / {m a +m b ), X a )S b (m b Pc / (m a +m b ), Xb). (12) 

This permits changing the integral over q into a delta function which sets r' to zero. Equation Ijlll) can then be 
written as 



dN^ 
d 3 Pn 



d 4 X a d 4 X b S a (m a P c /(m a + m b ), X a )S b {m b P c /{m a + m b ), X b )\u c {X' a - X' b )\ 



(13) 



This approximation is exact for thermal sources, since the Boltzmann factor is only a function of P a + P b , and is 
independent of the relative momentum. 

Unfortunately, the evolution matrix uc(X' a — X b ) involves a time difference. If the two times t' a and t' b are equal, u 
can be replaced with the relative wave function which is often a well-understood object. A third approximation can 
then be made by introducing time independent evolution matrix 0, i.e., 



\u c (X' a -X^^(27T) 3 \q>c(K-K 



(14) 



Without this last approximation it would have been essential to understand the details of the evolution of the particles 
for times between the emissions of the first and second particle. Since the particles are, by definition, not moving 
rapidly in this frame, the wave function at this time should not change substantially and this should be a reasonable 
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approximation. Putting the three approximations together, i.e., source factorization, smoothness, and ignoring of the 
offsets of emission times, we obtain an expression similar to that used in correlation studies 

dNi c ) r 

= (2tt) 3 d*X a d A X b S a {m a P c /{m a + m b ),X a )S b {m b Pc/{m a +m b ),X b )\<j> c (X' a -X b )\ 2 . (15) 



<PP, 



c 



After approximating the squared relative wave function as a 5 function, and noting that the phase space density 
can be related to the source function by the relation, 

f a (m a P c /(m a + m b ), x, t c ) = (2tt) 3 J d A X a S a {E ai m a P c /{m a + m b ), X a )6(x - X - v{t c - X )), (16) 

the coalescence formula of Eq. (0) would be reproduced, except that the source function in Eq. (|15J) is evaluated 
off-shell. 

m a E c /(m a + m b ) ^ E a (m a P c /(m a + m b )). (17) 

Within the smoothness approximation, there is no difference in replacing the energy arguments of the source functions 
in Eq. I|15|) with on-shell values. 

For a thermal source, the product of the source functions becomes a Boltzmann factor, exp(— Ec/T), and hence Eq. 
(|15fl is equivalent to a thermal expression. Thus, the coalescence expression, Eq. (0, the thermal expression, Eq. 10, 
and the source function expression, Eq. (|15|l . become equivalent for thermal distributions when the binding energy 
is zero. The lack of dependence on the binding energy in coalescence formulas is due to the sudden approximation, 
which presumes that the binding potential is turned on suddenly, i.e., energy is not conserved. From a quantum 
perspective, the source function in Eq. (|16|) is found by factoring quantum matrix elements with the density of states. 
If the matrix elements are assumed to be featureless, and if the density of states of the sources behaves as exp(_E/T), 
a thermal form is then obtained. When forming hadrons, it is hard to argue that the interactions with third bodies 
do not provide some additional leverage for creating final states with larger binding energies, i.e., the phase space 
density of protons should be greater than that of deltas. As can be garnered from the list of approximations above, 
neither the coalescence expression nor the thermal expression can be justified to better than the 10% level. Moreover, 
once a decision has been made regarding the choice of one of the formalisms, additional choices must be made for 
parameters such as the temperature and rest frame of the thermal source. For the coalescence formula it is necessary 
to make a choice of a rest frame for the coalescence. In this Lorentz frame the interaction between a and b turns 
on instantaneously. Momentum is conserved, but energy is not conserved. In all other frames, neither energy nor 
momentum are conserved. Again, it should be stressed that these subtle sensitivities to parameters and choice of 
reference frames disappear when the binding energy is small. 

To illustrate the sensitivity of coalescence pictures to binding energy, let us consider the extreme case of two 
constituent quarks of mass 350 MeV coalescing to form a pion. If the quarks are to coalesce in the two-particle rest 
frame, the energy will change by 60% from the coalescence in that frame. If the coalescence is observed in the lab 
frame, the pair will loose 60% of its momentum in the coalescence process. Whereas, if the coalescence frame is the 
laboratory frame, the laboratory momentum of the two quarks will be conserved in this process. Although, pions are 
an extreme example, most coalescence pictures involve mass changes at the 20% level, which provides a feel for the 
confidence with which one can apply coalescence prescriptions. 

Whether recombination is thermal or sudden by nature, the two expressions share an important property. Both 
formulas presume that the a and b components exist independently of one another. This will not be true if the 
partons, for example, originate from the same jet. It is the validity of the factorization that is related to whether 
or not hadronization can be described as the recombination of partons without regard to whether the two owe their 
existence to the same initial nucleon-nucleon collision. For many of the calculations presented in the following sections, 
the issue of thermal vs. sudden will be side-stepped by focusing on the coalescence of two quarks with constituent 
masses m q which combine to make a meson of mass 2m q . 

B. Coalescence and Elliptic Flow 

In this subsection, and in the examples presented in the next section, we consider the coalescence of n q quarks into 
a hadron of mass Mq = n q m q , where m q is the constituent quark mass. Thus, we ignore many of the subtle issues 
concerning the binding energy, and the relative merits of one formalism with another. The expression analogous to 
Eq. (|SJ) for the coalescence of n q quarks is 

f c (P, R, t) « [f q (m q P/M c , R, t)P . (18) 



5 



The hadron spectra is found by integrating the phase space density over coordinate space, 



d 3 P 



(25c + 1) 
(2 



^ J d 3 r[f q (P/n q ,R,t)] n \ 



(19) 



where Sc is the spin of the composite particle C. 

The principal goal of this paper is to investigate the behavior of the angular asymmetry which is quantified by the 
parameter v 2 defined in Eq. If the composite particle is made up of n q quarks, each with the same phase space 
density f q , we have from Eqs. Q and l|19|) 



,(™g) 



(Pt) 



/ d(j)d 3 r [f q (p t /n q ,ct),r,t c )} nq cos 20 
J d(/}d 3 r[f q (p t /n q ,(j},r,t c )} nq 



where <f> refers to the azimuthal angle of the momentum. 

The phase space density can be written as a product of the spectra and an effective density, p q (r, y t , i 
the relation, 



(20) 



defined by 



Pq{r,yt,<i>) 



fqiVt, (/>,*) 
(2Trf(dN(i)/d 3 p)' 



(21) 



Here, p q is non-zero over the region where particles are emitted. If the phase space density is constant within a fixed 
volume, SI, then p q = 1/S1. If p q is independent of <j), v 2 can be written in terms of the quark spectra, 



(n e ), , _ fdcf>[dNM/d 3 p] n «co S 2<f> 



Jd(j)[dN(iyd 3 p] nq 



(22) 



where dN^/d 3 p is evaluated for quarks at a fixed transverse rapidity, yt- If the quark spectra has an elliptic 
asymmetry as dN^ jd 3 p oc (1 + e cos 20), then v 2 q ^ — e/2 and w 2 can be expressed in terms of v 2 q \ 





1- 


f 2v ( 2 q) cos 24> 




cos 2(f) 


/# 


1 + 2v ( 2 q) cos 


2(b 


n q 



lEnLl,3,-("+l)(^ 



(n,-n)!([(n+l)/2]!) 2 



Zjn=0,2,- ^2 y (n,-n)![(n/2)!] : 



(23) 



For n q = 2, 3, the expression becomes 



,(2) 



2w 



(9) 



3w 



1 + 2(^)2 

( g) i+(4 g) ) 2 

l + 6(^ 9) ) 2 ' 



(24) 



For small w^'' one obtains the simple scaling relation, Eq. J2J. 

It is this simple scaling relation that has been cited as evidence for quark recombination However, the 

validity of this relation hinges on the assumption that the effective density is independent of <p. The nature of this 
assumption can be understood by considering two examples illustrated in Fig. ^ In the left side the effective volumes 
for emission upwards and to the right are identical, and the higher number of particles emitted to the right is due to 
the higher average phase space density of those particles. Thus, for this example, the elliptic asymmetry for dimers 
will be double the elliptic asymmetry for monomers and will satisfy the scaling relation, Eq. (J2J. For the emission 
illustrated in the right panel of Fig. ^ the average phase space density for emission in the two directions are identical, 
and the higher number of particles emitted to the right is due to a larger effective volume. In this case, the asymmetry 
for monomers and dimers is the same. 



FIG. 1: Two sources are illustrated above. The particles of a given velocity are represented by arrows. In (a) the effective 
volumes of the right-moving and upward-moving particles are the same but the right-moving particles have a higher phase 
space density. In (b) the phase space densities are the same but the effective volumes differ. In both cases there are twice as 
many particles moving to the right as upward so V2 is identical. But in (a) the V2 for a bound state of n q particles will increase 
by a factor of n q while in (b) the composite particle will have the same V2 as the constituents. 



To illustrate the same effect shown in Fig. Q] algebraically, let us consider a phase space density characterized by a 
Gaussian profile, 

f q (p t ,<fi,r) = 2 3 / 2 /(p t ,0)exp [~r 2 /(2R\p t A))} , (25) 

where / is the phase space density averaged over the coordinate space. The probability for emitting a coalesced 
n 9 -quark object is then 

dN c (2^ c + l)2^/ 2 n3/ n , 

W = (2^)3/* R <P"ti[f ■ (*» 

and t>2°' refer to the elliptic 

quarks 



If and v<£^ refer to the elliptic asymmetries of the average phase space density and effective volumes of the single 



(/) , * / d(j> f(p t ,4>) cos 20 



the asymmetry for the n 9 -quark object can be written as 

v^=vi n) +n q vi f) - (28) 

Equation l|28[) is again based on the assumption that the asymmetries are not large, but includes the chance that the 
asymmetries can originate from either asymmetries of the effective volume or asymmetries of the average phase space 
density. 

If the effective volumes are independent of 4>, is zero. The Eqs. <|23|) and (|28|) are then identical and the 
asymmetry scales linearly with n q . Such a case is illustrated in the left panel of Fig. ^ On the other hand, if the 

average phase space density is independent of 4>, v 2 is zero and the asymmetry is then independent of n q . This will 
be the case for emissions of the type illustrated in the right panel of Fig. ^ Several examples are explored in the 
next section, with the purpose of understanding what drives the asymmetry, the effective volume or the average phase 
space density. 

The effective volumes, R 3 in the Gaussian case, can be determined experimentally through correlation measurements 
l26l |27| or with coalescence measurements (2^, [2^, [2|j . This volume can be defined in terms of the phase space 
density, 

yt{y t ,4>) = r , 3 T- f , dj- 29) 

J d A r [f q {y t ,T) Y 
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For a Gaussian source of dimensions R x , R y and R z , fl — 2 3 / 2 R x R y R z . In considering the coalescence of nucleons, 
rather than that of quarks, the volume may be determined with two-nucleon correlation measurements. The effective 
volume may be also estimated from coalescence ratios. The coalescence volume can be defined by the ratio of spectra 
of two particles a and b that combine to form a coalesced object C as can be seen in Eq. (|29Jl . 

(2S C + 1) (dN a /d 3 p)(dN b /d 3 p) 
(2S a + l)(2S b + l) dN c /d3P ■ [ ' 

For small asymmetries, this yields ?4 = v 2 + v 2 ~ v 2^i which is basically a tautological restatement of Eq. 128fl . 
Other methods for measuring the source size as a function of the azimuthal angle, will provide independent means 
for determining whether the v^ 1 is independent of 0. Examples of such methods are identical pion correlations or 
proton-lambda correlations. 



III. BLAST WAVE MODELS 

Blast wave descriptions [lol ITU IT2I IT^ | , provide a simple means by which one can test thermal concepts along with 
those for collective flow. Typically, a blast-wave model requires four parameters: a temperature T, a breakup time 
r, a transverse radius R and a maximum transverse collective rapidity, pq. When considering elliptic asymmetries, 
one may also use two parameters for the transverse radii, R x and R y , and two parameters to describe the transverse 
collective motion. There are many variants of blast-wave models that differ from the form by which the space-time 
distributions are parameterized. Of the three parameterizations considered here, each assumes that the distribution 
of thermal sources is invariant to boosts along the beam directions. Although the choice of the form seems somewhat 
arbitrary when analyzing spectra, we find that the different forms result in rather different behavior in the elliptic 
flow. 



A. Shell parameterization for the blast-wave 



Asymmetries can be added to a thermal blast-wave in four independent ways via the temperature, the collective 
velocity, the chemical potential, or the volume. To compare the effects of these four types of asymmetries, we consider 
a simple azimuthally symmetric shell 



dN 
',(j) p dpt 



L dr] cosh7yp s (0„)exp 



- cosh 77 m tv /l + u]_{4>{s)) + u±((j> u )pt cos((j)p - 4> u ) 
T(4> u ) + nn(4> u )/T(4> u ) 



(31) 



Here u± is the transverse collective velocity, p(4> u ) is the angular density for sources having a collective flow in the 
radial direction <p u and 77 is the longitudinal rapidity of the source. The four asymmetries can be parameterized in 
the following way, 

p{4>u) = (1 + eycos2(/)„)p, 
exp[/a(0„)/T(0„)] = (1 + e^cos20 u )exp(n g /i/T), 
T(<f> u ) = (I + e T cos20 u )r, 
u±((f> u ) = (1 + e u cos2<j> u )u ± . (32) 

Inserting these equations into l|31|l . the expressions are then expanded to lowest power in e which allows the integrals 
to be performed analytically. The resulting expressions for each of the four expansion are: 



V2.V 



Jv 
Jo' 



V2„ 



V 2 ,T 



V 2 ,u 



Jo 




(33) 
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FIG. 2: Scaled values for elliptic flow as a function of scaled transverse momentum for constituent quarks each of mass m, = 350 
MeV coalescing to form composite objects of n q = 1,2,3 quarks with mass Mc = n q m q . Asymmetries are introduced into 
the cylindrical-shell parameterization of the blast-wave by adding elliptic distortions to the thickness of the envelope (upper 
panel), the density as driven by a chemical potential (second panel), the temperature (third panel) and the collective velocity 
(lower panel). Scaling is significantly violated for the upper two panels. 



where Jo, Je and J p are simple combinations of Bessel functions, 

Jo = 2K 1 (z E )Io(zp), 
Jv = K 1 (z E )h{z p ). 

Jp = -^K^Ze) [h{z p ) +I 3 (z p )] , 

Je = [K Q (z E ) + K 1 (z E )/z E ]I 2 (z p ), (34) 

with ze = (TOt/T)-\/l + u 2 ^ and z p — u±p t /T. 

Figure displays the resulting values of vi for each of the four asymmetries as a function of the scaled transverse 
momentum. For these calculations, the masses of the coalesced hadrons are assumed to be n q x 350 MeV. By plotting 
v-ijuq vs. the pt/n qi satisfaction of quark- number scaling of Eq. @ will provide indistinguishable lines for different 
values of n q . The upper panel shows asymmetries from adding an elliptic distortion to p{4> u )- This type of asymmetry 
will result from distorting the shape of the shell while keeping the magnitude of the collective velocity fixed, or by 
altering the thickness of the shell. Effectively, these distortions represent an elliptic asymmetry to the effective volumes 
of the source as described in Fig. ^ Scaling is then strongly violated since v 2 is approximately independent of n q . 
In the second panel, results are shown for adding an elliptic asymmetry to the chemical potential. This is same as in 
the upper panel, except that it scales with an extra power of n q since the density enters as exp(/z/T). Quark-number 
scaling then appears to be well satisfied except at low p t . 

The lower two panels of Fig. [21 show the result of elliptic distortions to the temperature T and to the magnitude 
of the collective velocity u. In both cases, quark-number scaling appears to be very well satisfied. But, in these 
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cases, satisfaction of quark- number scaling is driven by the fact that v 2 rises linearly with p t . If mass effects are 
neglected, v 2 can only depend on p t unless the chemical potential is varied, which will introduce an n q dependence. 
If v 2 rises linearly, scaling both the x and y axes by l/n q will yield a line with the same slope, and quark- number 
scaling should automatically be satisfied. In fact, universal curves will result if scaled by any number, but not just 
the quark number. It should be pointed out that the linear behavior with p t does not hold at low p tl where v 2 must 
be quadratic. The region over which it appears quadratic depends on the collective velocity u. For smaller u the 
quadratic region becomes noticeably larger. 

Even though the asymmetries in temperature and collective velocity produce excellent linear results for v 2 vs. p t , 
the superposition of several shells with different values of u and T can lead to non- linear behavior for v 2 vs. pt . This 
is because the slope can depend on the parameters, and by changing p t , the relative contribution of one shell may 
increase with respect to another. When v 2 vs. pt is not linear, quark-number scaling is more significant. 

B. Blast- Wave Parameterization of Retiere and Lisa 

The next parameterization we present was used to analyze elliptic flow and correlations in This model 

describes a uniform transverse density profile that is cut off by an elliptic surface. The normalized elliptical radius f 
corresponding to a surface of constant f is defined as 

f =/M- < 35 » 

Emission is confined to the region where f < 1 . The transverse rapidity is defined to rise linearly from the origin, 

Vt = f{p + Pa cos 20), (36) 

where the two parameters, po and p a parameterize the collective flow with p a driving the asymmetry. The transverse 
velocity is then u± = sinh(j/t). The direction of the collective flow is chosen to be perpendicular to the surfaces of 
constant f. 

The upper two panels of Fig. [3]show v 2 from this model with the parameters, T = 100 MeV, po = 0.9 and p a = 0.1. 
The calculations are for constituent quarks of mass m q = 350 MeV and with hadrons of mass n q m q . The results are 
independent of the proper time r and are not independently sensitive to R x and R y , but instead depend only on the 
ratio, R y /R x . If the matter falls apart quickly after the initial overlap of the colliding nuclei, then R y /R x > 1. The 
upper panel of Fig. [3] with R y /R x = 1.2 displays results for such a scenario. Here v 2 of the composite hadrons is 
calculated using the phase space density as described in Eq. (|19fl , and then v 2 is divided by the number of constituent 
quarks. The scaling criteria of Eq. J2J) is seen to be satisfied rather well in this case as all the curves lie close to one 
another. 

At later times, the unequal transverse expansion will overcome the initial out-of-plane extended profile and reach 
an in-plane-extended geometry. The second panel of Fig. [3] displays the results for R y /R x — 0.8. In this case, 
quark- number scaling is strongly violated. Moreover, v 2 becomes negative in some instances as was noted in |l3j. 
This peculiar behavior is due to the constraint that the collective flow is perpendicular to the constant f surfaces. 
Thus, non-zero v 2 values even occur in the case where p a = 0. 

The constraint that the collective flow is perpendicular to the surfaces of constant f is not justified from the 
perspective of hydrodynamics. It is expected that the accelerations, which are driven by the density gradients, are 
perpendicular to the surface of the ellipse, rather than the velocities. For early times, before the shape has a chance 
to react to the acceleration, the transverse acceleration and velocities should be parallel to one another. But at later 
times, when R y /R x approaches or exceeds unity, the constraint enforces a rather awkward velocity profile. This will 
be more apparent below where we consider the third parameterization. 

C. Hydrodynamically Inspired Blast Wave 

The second blast-wave model we will consider has the benefit of representing a solution to non-relativistic hydro- 
dynamics [12j |. This parameterization represents a solution to hydrodynamic equations in the limit that P = nT, 
and particle number is conserved. Hydrodynamic can realize this form even if the thermal motion is relativistic. The 
form also works for the simple ultra-relativistic equation of state, P oc e. However, this parameterization does not 
represent an exact solution in the limit of high collective velocities. 
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FIG. 3: Scaled values for elliptic flow as a function of scaled transverse momentum for constituent quarks each of mass 
m q = 350 MeV coalescing to form composite objects of n q — 1,2,3 quarks with mass Mc = n q m q . A constant temperature 
of T = 100 MeV is considered for all the calculations. The results are in the blast-wave model of Retiere and Lisa |13jl (upper 
two panels), and to a more hydrodynamically consistent parameterization (lower two panels). Quark- number scaling is better 
realized by the hydrodynamic parameterization. For the Retiere-Lisa parameterization, scaling is strongly violated when the 
ellipse becomes elongated in the in-plane direction. 



In this parameterization, the transverse profiles are Gaussian and the flow velocities rise linearly from the origin. 

N f ~ 2 - 2 



/ « 



(2nR 



-^-y exp^- - n 9 - j exp(-i?7T), 



(37) 



where u x and u y are the transverse relativistic velocities. As in the previous model, it is assumed that the matter 
disintegrates at a fixed temperature T. The factor n q in the exponential suggests that the space-time structure of 
the Gaussian is driven by a position-dependent chemical potential connected to the conserved quark number. In this 
parameterization the elliptic flow does not depend on the radii and collective flow parameters separately, but instead 
only on the combinations a x R x and a y R y . Since the collective velocities scale simply with the space-time coordinates, 
it is straight-forward to express the spectra over integrals of the collective velocities, 



dN 
d 3 p 



dr\ cosh j] u±du±d(j) u exp 



■ cos(0 u - <pp) 



2 

9 u x 



n q u 2 y 



T 



2a x R 2 x 2a v R? y 



(38) 



where <j) p refers to the direction of the momentum and <f> u is the azimuthal direction of the transverse collective velocity 
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Uj_. 

The integral can be performed analytically if the asymmetry is small. Defining the parameters u and e, 



1 

2 ( 

a x R x = u(l + e), ctyRy = u(l - e), (39) 



u = -{a x R x + a y Ry), 



leads to an expression for which the angular integrals can be performed analytically after expanding to lowest power 
in e. 



df] cosh 7/ u±du±d<j) u ( 1 + n g e "rj cos(20 u ) j exp — — yl + + — cos(0„ - <f> p ) 
The calculation of «2 involves an additional convolution with cos2<j} p , and gives 

v 2 = n q e J u ± du ± (ul/u 2 )e~ n ^/ 2a2 Kx{z E )h{z p ) / M , 
N = 2 f u ± du ± e- n ^ 2n2 K 1 (z E )I {z p ), 



(40) 



z E = (m t /T)^l + ul, z p = u ±Pt /T. (41) 

Results for this prescription are displayed in the third panel of Fig. |3 for the case of a 350 MeV constituent quark 
mass. The temperature T is considered the same as in the previous example and u is chosen to be 0.5. There is 
no peculiar negative elliptic flow from this parameterization, and constituent scaling appears to be maintained to a 
remarkably high level. Unlike the previous parameterization, we find no sensitivity to the ratio R y /R x . 

It can be shown that for small flow velocities and non-relativistic temperatures, quark-number scaling becomes 
exact. However, this limit is far from realized, and the apparent scaling is only approximate. 

To test the sensitivity of the scaling to the Gaussian shape of the profile, we repeated the calculation with an elliptic 
profile. For this calculation, the collective transverse rapidities are again chosen to rise linearly with the position, 
but the density is chosen to be uniform within an elliptic region. Also, as in the case of the Gaussian source, i; 2 
is independent of the spatial parameters. The results for the ellipse, shown in the lower panel of Fig. [3] are nearly 
identical to the Gaussian case. It thus appears that the crucial factor for scaling is that the collective velocity rises 
linearly with the position. 

The fact that both the elliptic and Gaussian profiles are fairly successful at reproducing quark-number scaling may 
be surprising given that the Gaussian profile for the collective velocities has an explicit dependence on n q , while the 
elliptic profile does not. Since the Gaussian profile behaves as exp(— n q u 2 /2u 2 ), the average collective flow decreases 
as Lower collective flows result in lower values for i?2. However, since the distortions scale linearly with the 

density, ~ (1 + n q ecos2<fr u ), the two effects somewhat cancel one another to some extent, and one obtains values of 
V2 which depend mainly on p t . Since this dependence is roughly linear, scaling both the x and y axes by n q , or any 
other constant, still results in a semi-universal curve. 

Hydrodynamics provides for many more possibilities than the Gaussian density profile with linear velocity profiles 
as described above. A strong first-order phase transition will result in shock waves [H E3, which can strongly 
violate the assumption of linear velocity profiles. Additionally, the temperatures might not be independent of the 
position. Nonetheless, most full hydrodynamic solutions do tend to approach semi-linear forms by the time breakup 
is approached [TH . 

D. Sensitivity to Mass 

In addition to the assumption that the effective volume referred to in Eq. (|28() was independent of <fi, the scaling of 
i>2 with quark number requires that for a particle C composed of particles a and 6, 

fc{M c u, x) = f a (m a u, x)f b (m b u, x), (42) 

where u = 7V is the relativistic velocity. For a thermal source where 

fi(p,x) = vcp[-E'jT{x) + (jh(x)/T(x)], (43) 

Eq. I|42|l is satisfied with a thermal distribution if the three criteria are satisfied: 
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1. The temperatures describing a, b and C are the same. 

2. The energies are conserved, Ec(Mc) ~ E a {m a Pc / (m a + mb)) + Eb{mbPc / {m a + mb)). 

3. The chemical potential for C is the sum of the constituent chemical potentials, fic(x)/T(x) — fi a (x)/T(x) + 
fj, b (x)/T(x). 

The first criteria requires only the assumption that the different species are kinetically thermalized. 

Since the energies are a mass multiplied by the gamma factor describing the motion of the particles in the frame 
of the source at x, the second criteria can be re-stated as a requirement that the masses sum, Mc = m a + wj. This 
is equivalent to a requirement that there is no binding energy involved. However, even if binding energy is involved, 
another simple scaling can be derived, as discussed by Molnar . If the chemical potential is ignored, Eq. l|4^jl may 
be written as 

/ M =o(m a u,x) = [/ fI=0 (mou, I )] mi/m », (44) 
which then implies that if one can ignore the contribution to V2 from the volume like term, 

4 C) = ^ Q) = ^ b) . (45) 
m a m b 

For the case where the quark masses are the same, this results in 

4 C) (V) = ^ 9) (V). (46) 



rn 



The scaling with mass differs from the quark-number scaling in two ways as can be seen by considering the case of 
pions and protons. In quark number scaling the proton will have 50% more vi than a pion, while the difference will 
be approximately a factor of seven for the case of mass-number scaling. Secondly, for coalescence arguments based on 
conservation of momentum in the lab frame, various particles are to be compared at the same value of p/n q , whereas 
for mass-number scaling, the same velocity is to be compared, i.e., at the same value of p/m. Thus, for the case of 
mass-number scaling, a universal curve will be obtained by plotting Vijm vs. pt/m. Experimental results seem to 
more closely follow quark-number scaling. 

In order to understand the degree to which mass effect destroy the scaling of V2 with quark-number, we repeated 
the calculations of the previous subsection for the hydrodynamically inspired blast-wave model. Instead of using 
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hadronic masses equal to M c = 350n q MeV, we use the pion mass for the two quark state and the proton mass for the 
three quark state. The comparison of the two choices for masses is displayed in Fig. Remarkably little variation 
is observed for the different curves. We also found no observable changes for the shell-profile blast-wave calculations 
with similar mass modifications (not shown here). A significant sensitivity to the mass is however observed if one 
focuses at low p t below 500 MeV/c. The sensitivity is also more apparent for lower collective velocities. 

The third criteria for implementation of the thermal expression, i.e. the addition of chemical potentials, implicitly 
requires that there are conserved charges or numbers with Qc = Q a + Qb- If quark number is conserved in the 
hadronization process, the chemical potential can be associated with the net number of quarks and anti-quarks. This 
is analogous to the implicit assumption in coalescence that the quark content is unaltered by hadronization. Whereas 
this assumption is certainly reasonable for the coalescence of protons into deuterons, it can be questionable for the 
case of quark coalescence since quark anti-quark pairs can be copiously created in hadronization as expected either 
from entropy arguments or from dynamic pictures that involve string breaking. However, for high p t mesons, it does 
seem reasonable to neglect fragmentation processes, since such processes involve a color flux tube to form between 
the outgoing quark and the thermal medium. This flux tube lowers the momentum of the outgoing quark by roughly 
one unit, and in an exponential spectrum, can be neglected at values of pt well above the temperature. The breaking 
of such strings should be important for observables concerning particles with p t < 1.0 GeV/c. One can make strong 
arguments that hadronization for quarks below 1 GeV is of much different character than the hadronization at higher 
values |3|. 

IV. MICROSCOPIC SIMULATIONS 

More sophisticated models of heavy ion collisions, such as hydrodynamic or Boltzmann descriptions, provide more 
complicated forms for the final phase space density. By analyzing the results of these numerical models, one can 
gain insight into the validity of the assumptions of blast-wave parameterizations. Unfortunately, there is not yet a 
"standard model" for the evolution of a RHIC collision. Hydrodynamic models using equations-of-state based on 
lattice gauge theory tend to significantly over- predict the lifetime of the collision [2^ |2y, • Boltzmann approaches 
which effectively represent stiff equations of state tend to somewhat under-predict source sizes |34j . Since the issue of 
the space-time development of the collision remains uncertain at the 50% level, we will instead focus on a simplified 
Boltzmann description. By adjusting the cross sections and the initial conditions, we intend to better understand 
how and why ^2(2/*) can depend on n q . 

For the first calculation, we consider a system of 500 constituent quarks per unit rapidity of mass 350 MeV and 
initially thermalized at a temperature of 275 MeV. The particles are then allowed to collide with simple s-wave 
cross sections using the code GROMIT |35J to solve the Boltzmann equation. The initial density distribution in the 
transverse directions are set by either a Gaussian distribution, or by an ellipse with constant density. The initial 
sizes for the Gaussian profiles are chosen to be R x — 2 fm and R y = 2.6 fm, while for the ellipse, the radii are 4 and 
5.2 fm. The longitudinal distributions are initialized according to boost invariant criteria, v z = z/t. Boost-invariant 
constraints are maintained by cyclic boundary conditions along the z axis. 

The upper panel of Fig. [SJshows calculated values of V2 from simulations using three cross sections, 5 mb, 20 mb and 
80 mb. In addition to the vi for constituent quarks (n q = 1), V2 for coalesced mesons (n q — 2) are also shown. Meson 
coalescence is calculated by sorting the quarks into bins according to their transverse momentum p t and azimuthal 
angle (j) p , then coalescing them with a weight proportional to exp(— 5r 2 /2r 2 ), where Sr is the separation of the two 
quarks as measured in the two-particle rest frame. The coalescence distance r c is chosen to be 1 fm. Constituent 
number scaling appears to be satisfied to high precision for all three cross sections. This is probably due to the fact 
that the initial conditions and equation of state satisfy the conditions for the scaling hydrodynamic solution discussed 
in the previous section, which in the limit of small flow velocities can be shown to lead to quark-number scaling. 

Calculations based on an initial elliptic profile with sharp cutoff are displayed in the lower panel of Fig. [5] Aside 
from the density profile, the calculation was the same as was shown for the Gaussian profile with a cross section of 
20 mb. In this case, quark-number scaling is significantly violated. This failure can be attributed to the shock wave 
which develops from the discontinuity in the density. From the perspective of hydrodynamics, there is no pressure 
gradient except at the surface. This forces the development of a shock wave. Emission is then more surface-like and 
less volume-like, which is more suggestive of the illustration on the right hand side of Fig. ^ F° r surface-like emission, 
one expects the asymmetry to owe itself to the disparity in the effective surface areas for emitting in- and out-of-plane, 
rather than to asymmetries in the respective phase space densities. For equations of state that incorporate a large 
latent heat, shock waves also develop since there is a large region of mixed phase for which there is no pressure 
gradient and therefore no transverse acceleration |30l l3l| . One may also expect to observe violations from scaling in 
these instances. It will be interesting to analyze hydrodynamic calculations to see whether the presence of a large 
latent heat will, indeed, manifest itself in the behavior of v-i- 
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FIG. 5: Scaled values for elliptic flow V2/n q as a function of scaled transverse momentum pt/n q from Boltzmann transport 
calculation. Upper panel: The initial state evolved from Gaussian density profile with elastic scattering cross sections of 
a — 5, 20 and 80 mb (bottom to up). Nearly perfect scaling with quark number is exhibited. Lower panel: The initial state 
evolved from elliptic density profile with a sharp cutoff. V2 significantly violate quark-number scaling. 



V. SUMMARY 

Quark-number scaling for elliptic flow has been advocated as a signal for quark recombination. The two principal 
goals of this study are to identify precise conditions for quark number scaling, and to determine the degree to which 
thermal, coalescence, and microscopic models may realize these conditions. 

Two conditions are necessary for realizing quark-number scaling. First, the phase space density for a composite 
object must equal the product of the phase space densities of the constituent quarks. Second, the effective volume for 
quarks that have a given momentum must be independent of azimuthal direction. It will be much more difficult to 
determine the validity of the latter condition since the effective volume is determined both by details of the spatial 
source size and by the collective flow and temperature. 

The requirement of factorization is manifestly satisfied by coalescence models. It is approximately satisfied by 
thermal models, and becomes exactly satisfied in the limit that the masses of the coalesced hadron equal the product 
of n q and the constituent quark mass. However, we never observe a significant dependence to the behavior of V2 as we 
varied the masses of the final hadrons for pt/n q I 1 GeV/c, the range of interest for quark recombination arguments. 
It should be emphasized, that in a thermal model, the only means to have non-uniform densities are by varying 
the temperature and chemical potential. If non-uniform chemical potentials are responsible for the density profile, 
thermal models implicitly assume that charge must be conserved by hadronization. For the hydrodynamically inspired 
blast-wave with a Gaussian profile, the temperature is kept constant, and the density is driven by non-uniformities in 
the chemical potential. Quark- number scaling for this model is then dependent on the baryons having 3/2 the charge 
of a meson, although this sensitivity is weak due to the linearity of V2 vs. p t . 

The models invoked in this study are schematic by nature, but they provide several important informations regarding 
the interplay of coalescence and elliptic flow. It does appear that linear velocity profiles with sudden time- like break-up 
surfaces promote constituent-number scaling criteria of Eq. @ ■ This is verified by both blast- wave and microscopic 
modeling. Since linear velocity profiles tend to develop, except in the presence of strong shocks, Eq. should be 
valid unless the initial density profile is abnormally sharp, or if the equation of state has a strong first-order phase 
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transition. 

Furthermore, the scaling condition can be tested by measuring the source-size as a function of the azimuthal 
volume. This can be done either with correlations or with coalescence ratios. If the source-size volumes at a given 
pt are independent of the azimuthal angle, it follows that constituent-number scaling of vi should be satisfied if the 
factorization criteria are met. Perhaps the easiest way to test this will be to estimate the d/p 2 ratio for deuterons 
and protons at the same transverse rapidity. This ratio is inversely proportional to the source size volume, so if 
constituent-number scaling is valid, it should be independent of the azimuthal angle. In other words, the vi for 
deuterons will then be twice v-i for protons. 

We have explored numerous models for the satisfaction of quark-number scaling: 

1. Four different distortions of a blast- wave with an elliptic-shell geometry. 

2. The blast-wave model of Retiere and Lisa, with both in-plane and out-of-plane extended profiles. 

3. A blast-wave with simple linear behavior for the collective velocity and either a Gaussian or sharp elliptic profiles 
for the density. 

4. A solution of the Boltzmann transport equation with a coalescence prescription for combining quarks into 
hadrons. Both Gaussian and sharp elliptic profiles are used for the initial densities. 

In each of these prescriptions, the factorization condition is explicitly satisfied. Thus, the satisfaction of quark-number 
scaling depend solely on whether the effective volumes are ^-independent. Unfortunately, it is not algebraically possible 
to show that any of the models explicitly satisfy quark-number scaling, except for the blast-wave with linear velocity 
gradients and a Gaussian profile in the limit of small expansion velocities and non-relativistic temperatures. However, 
an empirical pattern does seems to emerge: Linear velocity gradients and sudden volume-like disintegrations seem to 
closely satisfy quark-number scaling. Hydrodynamic evolutions tend to approach this geometry, except in the case 
of surface-like emission, which will be favored if shock waves are to develop. Shock waves are expected if the initial 
density has a sharp fall-off, or if there is a large latent heat which will lower the interior pressure of the fireball to 
give a shock wave the opportunity to develop. 

Quark recombination is synonymous with satisfaction of the first condition, factorization of the phase space density 
of the combined object. This is because factorization assumes that two-particle densities are products of single-particle 
densities, i.e., jet-like correlations are lost. However, since factorization represents only one of the two conditions for 
quark- number scaling of i>2, there is always some ambiguity in concluding that both conditions are satisfied when 
one observes the scaling. For instance, it could be that neither condition is satisfied, but that the failures to meet 
both conditions conspired in some way to mimic the satisfaction of both criteria. Of course, a stronger case can be 
made if scaling is observed for a large number of cases. However, one must also remember that if the various cases 
involve hadrons with different cross sections, i.e., strange vs. non-strange hadrons, then one may expect violations of 
V2- This can be understood by viewing Fig. [3] where the overall value of V2 is shown to be strongly sensitive to the 
cross section. 

Our investigations of both model results and of the underlying theory shows that observing quark-number scaling 
of i>2 indeed suggests that quark recombination may be at hand, but there are quite a large number of caveats that 
prevent one from declaring this as a stand-alone signature. However, it is clear that V2 is, at the very least, very 
sensitive to the evolution, break-up mechanism and hadronization of the fire-ball. Thus, the details of the behavior of 
V2 can be used to disqualify a large number of potential models, even if pointing to the root cause of the disqualifying 
attributes might be somewhat ambiguous. 

We conclude by recommending that source-size measurements at intermediate pt can significantly clarify some of the 
issues surrounding quark recombination. First, interferometric or coalescence measurements may determine the source 
size as a function of pt and </>. If source sizes are seen to be independent of the azimuthal angle, it should support the 
statement that quark-number scaling demonstrates that phase space densities factorize. More importantly, a simpler 
and more direct point can be made about such source-size measurements. If all particles in the intermediate-pt region 
with the same azimuthal angle originate from the same jet, the source size from correlation or coalescence analyses 
will be identical to that expected from pp collisions of ~ 1 fm. If the emissions are independent, and factorization is 
valid, the source sizes will reflect the entire emitting region of such particles. Even if emission of high energy particles 
is confined to the surface during the first few fm/c, the effective volumes will be on the order of tens of cubic fm. 
The real value probably lies somewhere in between. Since jet- like correlations are observed at these pt region |36j| . 
the assumption of completely independent sources is probably not viable and more detailed models which incorporate 
jet-like correlations into coalescence pictures are currently being developed H3,|3£|. A measurement of the source 
volume will then provide crucial quantitative information to decide the relative role of recombination. 
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